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Abstract. 

We study the stabiHty properties of isolated star forming dwarf galaxies 
which undergo dynamically driven star bursts induced by stellar feedback. Here 
we focus on the impact of the adopted ISM model, i.e. either a diffuse or a 
clumpy ISM. We apply a one-zone model extended for active dynamical evo- 
lution. This allows for a coupling between the dynamical state of the galaxy 
and its internal properties like star formation activity or the thermal state (or 
dynamical pressure, respectively) of the interstellar medium (ISM). 

We found two major types of repetitive star bursts: one set (type A) of 
quasi-periodic star bursts is related to the dynamical timescale of the galaxy. In 
that case, the star formation follows the variations of the gas density induced 
by decaying virial oscillations. The second set (type B) of recurrent star bursts 
is characterized by a long quiescence period given by the sum of the dynamical 
and the dissiptional timescale: after a first burst, the inserted energy leads to a 
substantial expansion of the system, by this stopping any significant SF activity. 
A next burst might occur, when the gas reaches high densities again, i.e. after 
the gas recollapsed and the energy injected by stellar feedback is dissipated. 

In case of a diffuse ISM model, type A bursts are the most common type 
due to the high efficiency of radiative cooHng (no type B bursts are found). 
Bursts occur then mainly during an initial transitory phase. In case of a clumpy 
ISM model (i.e. dissipation by inelastic cloud-cloud collisions), the dissipative 
timescale is of the order of the dynamical time or longer. This allows for both, 
type A and type B bursts. Whereas initial transitory bursts are quite common, 
type B bursts are only found in a small mass range for given feedback and 
dissipation parameters. 



1. Introduction 

Starb ursts are a well known phenomenon in some dirrs (e.g. Gallagher & Hunter 
( 19841 ): for a critical discussion of the exact definition of starbursts see Knapen & 
James ( 20091 )). They imply a mechanism to organise large-scale star formation 
which cannot operate all the time due to fast gas consumption. Moreover, the 
instability leading to an episodic SFR should work in a small parameter range 
only, because most of the isolated dIrrs seem to evolve in a highly self-regulated 
manner. 
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The first theoretical models dealing with the variability of the star formation 
rate in galaxies were based on closed-box models. For reasonable models of the 
interstellar medium (ISM) it turned out that stell ar fee dback is rather efficient in 
suppr essing instabilities (e.g. Ikeuchi &; Tomita ( 19831 ) . Scalo & Struck-Marcell 
( 19861 )). These models only allow for a burst-like behaviour when a long time- 
del ay bet ween star formation and stellar feedback is introduced. Koppen et 
al. |l995l . hereafter KTH95) corroborated this result by investigating a simple 
numerical model motivated by full chemodynamical simulations. They demon- 
strated that the effective SFR is almost independent of the detailed recipe for 
the stellar birth function, provided a negative feedback due to the thermal state 
of the ISM is considered in the stellar birth function. 

These box models, however, suffer from neglecting galactic dynamics. Sev- 
eral attempts have been made to study the evolution of isolated dwarf galaxies 
by taking full dynamic s, as well as star formation and s tellar feedback, into ac- 
count (Kruegel et al. ( 19831 ) . and later Hensler et al. ( 20041 ). Pelupessy et al. 
(|2004l 1. Struck (|2005l 'l). Though these calculations differ in many details, they 
agree in producing large-scale star formation variations on dynamical timescales 
which indicates the crucial role dynamics plays for episodic starbursts. 

A major disadvantage of numerical simulations is their complexity. This of- 
ten prevents a detailed investigation of the parameter space. Instead of analysing 
dwarf galaxies by such detailed models, we apply a set of equations that is as 
simple as possible, but still incorporating the generic features of the complex 
numerical simulations. We use an extended version of the analytical models 
of KTII95 that allows us to deal with the dynamical evolution of a galaxy in a 
simple way (for details cf. to Theis & Koppen (20091 )). Here we compare the evo- 
lution of dwarf galaxies for two different ISM models: in the diffuse ISM model 
energy is dissipated by radiation, whereas in the clumpy ISM model energy is 
lost via inelastic cloud collisions. 



2. Results 

Diffuse ISM models. Fig. [1] shows the evolution for the diffuse ISM model. 
Due to strong cooling, the gas temperature drops almost immediately to values 
near lO'^K bringing the system far out of virial equilibrium. In the ongoing 
collapse the density increases and the gas temperature drops even more. The 
cooling and density enhancement continue, until star formation becomes more 
prominent after about 100 Myr. The first stars reinject energy to the ISM 
leading to the quasi-equilibrium stage known from the box models of KTII95. 
However, different to them, the density still increases due to the ongoing collapse. 
Though the star formation and, by this, the stellar energy injection grows, the 
feedback cannot stop further collapse, because cooling is too efficient to prevent 
the gas from reaching virial equilibrium. The collapse proceeds until the angular 
momentum barrier becomes dominant and the system expands. Then the first 
star burst episode finishes and the gas is quickly heated to lO^K. However, due 
to the strong enhancement of the cooling rate at temperatures beyond lO^K, 
the gas never reaches the virial temperature. 

The longterm evolution of this model shows that the previously discussed 
bursts are just transitory phenomena. The oscillations are almost completely 
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damped after 2 Gyr. Such a behaviour, i.e. only an initial burst, was found in 
almost all models with radiative cooling as source of dissipation. 
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Figure 1. Evolution of a model with a baryonic mass of 2 ■ 10^ M©. Shown 
are the masses of the baryonic components (upper left), the mean radius of the 
system (upper right), the total SFR (lower left), and the gas temperatures, 
i.e. the actual temperature and the virial temperature (lower right). 
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Figure 2. Temporal evolution of a model with energy dissipation by inelastic 
cloud collisions. The baryonic mass amount to 4.5 • 10^ M0 (for details see 
Fig. [ID. 

Clumpy ISM models. The evolutionary behaviour changes when applying the 
clumpy ISM model. Though the structure of the equations remains identical, 
the dissipational timescale differs substantially from the previous case in which 
the cooling timescale was much shorter than other involved timescales. For 
cloud-cloud collisions, the timescale Tcou becomes of the order of Gyrs. 
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Driven by stellar feedback and the much longer dissipational time the ISM 
expands after the first burst (Fig. [2l upper right). The dark matter halo, how- 
ever, prevents the gas from being blown-away. Therefore, it collapses again 
reaching a dynamical equilibrium after about 2 Gyrs. This can be seen from a 
comparison between the actual kinetic temperature (velocity dispersion of the 
cloud system) and the virial temperature of the system. The ongoing evolution 
is governed by a slow loss of energy due to dissipation. During this stage of 
slow shrinkage the evolution is close to virial equilibrium. After 6 Gyrs the sys- 
tem lost sufficient energy to start with the next run-away collapse resulting in 
another starburst event after a long period of quiescence. Finally, the galaxy 
produces stars on a mean rate of about O.lMQyr"^. 

Discussion 

The main difference between the two model scenarios is the ratio r between the 
dissipational timescale t^iss and the dynamical timescale r^yn- In case of the dif- 
fuse ISM model (small r) the energy input by the stars is almost instantaneously 
radiated away. Therefore, the radial oscillations are only decaying initial virial 
oscillations acting on a dynamical timescale. This behaviour turned out to be 
rather robust against variations of model parameters like stellar heating, gas 
mass, initial size or stellar birth function. 

On the other hand, when r becomes of the order of one or larger (i.e. for 
the clumpy ISM model), two timescales are important: the dynamical timescale 
and the dissipational timescale. After a first burst, a long period (rdyn + Tdiss) 
of quiescence follows, before the next star formation episode starts. 

During most of the quiescence phase the gas is semi-adiabatically accreted 
with an average infall velocity of only a few kms~^. The related extended gas 
reservoirs are known around many dwarf galaxies (e.g. NGC 4449). They are 
very sensitive to perturbations by other galaxies. E.g. a substantial amount of 
the gas might be stripped off, by this reducing the strength of the second burst. 

After the second burst the galaxy enters a phase with continuous star forma- 
tion up to the present (here: 10 Gyr). This is in agreement with observations of 
several dirrs which still show a substantial, more or less con tinou s star formation 
ongoing over the last one or few Gyrs (e.g. Greggio et al. ( 19931 )). 
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